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Abstract 

Using two specific models and a model-independent formalism, we show that 

in addition to the usual quadratic "side", "out" and "longitudinal" terms, 

a previously neglected "out-longitudinal" cross term arises naturally in the 

exponent of the two-particle correlator. Since its effects can be easily observed, 

such a term should be included in any experimental fits to correlation data. 

We also suggest a method of organizing correlation data using rapidity rather 

than longitudinal momentum differences since in the former every relevant 

quantity is longitudinally boost invariant. 
25.70.Pq 
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Typeset using REVT^X 



The experimentally measured Hanbury-Brown Twiss (HBT) correlation between two 
identical particles emitted in a high energy collision defines a six dimensional function of the 
momenta pi and P2 [Q- A popular way of presenting these is in terms "size parameters" 
derived from a gaussian fit to the data of the form 



C(q, K) = 1 ± A exp [-q 2 s R 2 s (^) ~ ?X(K) - <?^ 2 (K)j (1) 

where q = Pi — P2, K = ^(Pi + P2), the + (— ) sign is for bosons (fermions), and the HBT 
cartesian coordinate system is defined as follows ||: The "longitudinal" or z (subscript I) 
direction is parallel to the beam; the "out" or x (subscript 6) direction is parallel to the 
component of K which is perpendicular to the beam; and the "side" or y (subscript s) 
direction is the remaining direction. 

In this letter we assert that significantly more can be learned and better fits achieved if 
an "out-longitudinal" cross term is included in any gaussian fits to the data. In other words, 
we suggest that the data should be fit to a function with the following form 

C(q, K) = 1 ± A exp [-q 2 s R 2 s (K) - q 2 R 2 (K) - q 2 R 2 (K) - 2q ogi R 2 ol (K)] , (2) 

where R? ol is a parameter which can be either positive or negative; we simply use the R? 
notation to denote the fact that it is has the dimension of an area. 

Since particles 1 and 2 are indistinguishable, the overall sign of q is irrelevant. The 
relative signs of the various components of q, however, are well defined physical quantities 
for any given pair. Our sign convention for q will be such that q s is always positive. We 
can thus unambiguously discuss correlations for negative as well as positive values of both 
q Q and q t . 

To see how an "out- longitudinal" cross term arises in two-particle correlations, we use 
the following well established theoretical approximation 

C(q,K)^l± Ud , xS{XiK)? (3) 



where q^ = Ei — E 2 , K = Ek = \Jm? + |K| 2 . Here S(x,K) is a function which describes 
the phase space density of the emitting source. For pairs with |q| <C Ek 



q-x ~ (P q + Piqi) t - q x - q s y - q t z , (4) 

where $ = Ki/E K . 

As a simple example, we consider the following cylindrically symmetric gaussian emission 
function 



S{x,K) = f{K)e W 



x 2 + y 2 z 2 (t-t ) 2 



(5) 



2R 2 2L 2 2(6t) 2 

Using (§) and ([|), it is easy to see that the corresponding correlation function takes the 
form: 

C(q, K) = 1 ± exp [-q 2 R 2 - q 2 Q (R 2 + (3 2 (6t) 2 ) - q 2 {L 2 + (3?(5t) 2 ) - 2q q l p o p l (St) 2 } . (6) 

For this model, the q Q qi cross term thus provides a measurement of the duration of particle 
emission (St). 

One might think that this cross term is just a trivial kinematic effect which would not 
arise if the correlation were calculated in some more carefully chosen coordinate system or 
reference frame. For example, for spherically symmetric systems (L = R), it has been shown 
that the cross term vanishes if axes are chosen parallel and perpendicular to K (rather than 
parallel and perpendicular to the beam) || . The reader can verify, however, that for systems 
with L ^ R, the cross term does not vanish in these rotated coordinates, but rather measures 
the L 2 — R 2 asymmetry of the source. 

Another system that is often proposed is the Longitudinally Co- Moving System (LCMS), 



which is defined as the frame in which /?/ = |I|f|,|5|. Glancing at (|6|), it naively appears 
that the cross term will vanish in this frame. Being more careful, however, one can see that 
after transforming to the LCMS frame (primed variables) 

t = ji(t-/3 l z) z' = ^(z - ftt) , (7) 



where 7/ = l/Jl — (3f, t'z' terms arise in the transformed emission function 5". These in 
turn lead to a nonvanishing cross term of the form 
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Kx = /W [(Sty + l 2 \ , (8) 

where Pi and 7; are evaluated in the CM frame. We can see that the cross term cannot 
in general be removed simply by transforming to another coordinate system or reference 
frame, and therefore it is certainly not just a trivial kinematic effect. In fact, R 2 { contains 
physical information about the emitting source which is just as important as that found 
from evaluating the difference R 2 Q — R 2 S . 

In order to get a broader feeling for what the cross term measures, we now introduce 
a general formalism valid for any cylindrically symmetric emission function which can be 
expressed in a roughly gaussian form. Using @, we expand exp(iq-x) in (§) for q-x < 1 to 
find (see Ol) 



C(q, K) = 1 ± {l - q 2 (y 2 ) - ([q (J3 t - x) + q x {fct - z)f 

+ (q Q (J3 t -x)+ qi (At - z)) 2 + O [(q-x) 4 ] | , (9) 

where the q s q Q and q s qi terms vanish due to cylindrical symmetry || and we have introduced 
the notation 

<« - <«<k) - s f£li:% • ^ 

Exponentiating (Q), we can see that for any cylindrically symmetric system, the corre- 
lation function for small momentum differences (qiRi 1) can be expressed in the form of 
Eq. (0) where A = 1 and the functions R 2 (K) can simply be read off as the coefficients of 
the corresponding qiqj terms in Eq. (|]): 



R 2 — 


(y 2 ] 


) - (y) 2 = 






Rl = 


((x 


- Pot) 2 ) - 


(X - Pot) 2 




R 2 = 


«* 


- Pit) 2 ) - 


(z - Pit) 2 




R li = 


((x 


-p o t)(z- 


Pit)) - (X- P Q t)(z- Pit) . 


(11) 



These expressions have the typical form of variances and show that the HBT "size param- 
eters" are really lenghts of homogeneity associated with the source function S(x,K) ||12|| . 



For a static source, these homogeneity lengths are equal to its geometric size in the various 
directions which can then be directly extracted from the HBT correlator. For expanding 
sources the interpretation of the HBT size parameters is more involved |f?| ^T0] , [T~2| , [T3]1 , and 
the HBT parameters are usually smaller than the geometric extensions of the source. The 
cross term is seen to measure the temporal extent of the source as well as the x z, xt, and 
z t correlations of the emission function. Note that the LCMS radii can be found from (|TTD 
by setting (3\ = and using S(x',K') (see Eq. (0)). The cross term vanishes in this frame 
if and only if the source S(x' , K') is reflection symmetric under z' — > —z' (which is not the 
case for our source (§D). 

One might argue that the model independent expressions of Eq. (|XI|) should not be 
compared to experimental correlation radii since the former measure second derivatives of 
the correlation function around q = (because we used q-x C 1 to derive them), while 
the latter are parameters of a gaussian fit to the whole correlation function @-|5|]. On the 
other hand, for any source which has a roughly gaussian profile in some complete set of 
spatial coordinates, the two different methods of measuring radii will give roughly the same 
results. For these types of models, the simple expressions generated by Eq. ([□]) provide 
valuable insights as to how various parameters of the source distribution qualitatively affect 
measurable features of the correlation function. 

We already discussed one gaussian model in Eq. (|5|), but here we would like to discuss 



another, possibly more realistic model which is similar to the ones presented in ||13|| . In the 
center of mass frame of an expanding fireball, we define the following emission function 

K-u(x) p 2 rf (r — r ) 2 



S(x,K) = . exp 

(27r) 3 r v /27r( ( 5r) 2 



T 2R 2 G 2(5r]) 2 2{8rf 



(12) 



where T is a constant freeze-out temperature, p = ^/x 2 + y 2 , r = y/t 2 — z 2 , 
i] = | ln[(t + 2;) / (t—z)), m t = \jm 2 + Kj_, and Y is the rapidity of a particle with momentum 
K. Note that in the limit 5t — > 0, (0) becomes the Boltzmann approximation for a 
hydrodynamic system with local flow velocity u(x) which freezes out on a 3-dimensional 



hypersurface of constant longitudinal proper time Tq and temperature T ||14|,p|. We will 



consider a flow which is non-relativistic transversally but which exhibits Bjorken expansion 
longitudinally 



u(x) - ( (l + \{vp/R G f) ch v , (vx/R G ), (vy/Ra), (l + \{vp/R G ) 2 ) sh v 



(13) 



where v <C 1 is the transverse flow velocity of the fluid at p = R G . In [BJ we show that 
when the emission function (^) is integrated over spacetime, it produces a very reasonable 
one-particle distribution. 

To facilitate calculating the correlation function, we can make the physically reasonable 
assumption that 5t/tq < \ so that we can be justified in replacing integrals over only positive 
values of r with ones ranging from — oo to +00. We can then achieve analytic results by 
making a Taylor expansion || in the parameter 
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1 m t 



(14) 



(Sv)l (Sr)) 2 T ' 

Note that for pairs in which m t /T ^> l/(<5?7) 2 as were studied in |TJ|], (5rj) 2 becomes simply 
T/m t . 

Using the expressions ( |TT| ) and keeping a few subleading corrections which are important 
when considering pions [El, we find 



Ri 



R~ 



Rl + 



K 2 K 2 K 2 



mi 



-(Sry + ^(3fT 2 (6 V )i + 



mi 



+ 



K 2 



rn 



2 T o 



(3 2 v - 2A 



{5 V ) 2 2 



mi 



(Sv) 



Y 



(Sr) 2 (5r,)l 



4 . 

* ) 



^r 2 (5 V )l + ^t(5t) 2 (5 V )1 + ^r 2 (5 V )t 



El 



E 2 K 



E 2 K 



R. 



ol 



-M3iT 2 {8rj)l-p o 



0i 



Y 



(5r) 2 (5r ] )l - ftr, 



Y 



where 



(15) 



1 1 / m t 2 
Rj = R%{ l + Y V 



(16) 
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and v = 1 + (R*/Rg) 2 — k{ m t/T)(Sr])l. As pointed out in |L3| and seen from Eq. ( |Tt^ 



transverse flow causes the "side" radius to measure something smaller than the real geomet- 
rical radius Rq- What it does measure is the transverse region of homogeneity of the fluid 
as seen by particles with a given p t It is also interesting that R 2 a — R 2 S depends on the 

average rapidity and is not quite directly proportional to the duration of particle emission 
St even for pairs with f3i = ||16|| . In our opinion, however, the most interesting feature of 
this model is the cross term radius R 2 ol . Although just as in Eq. @ the cross term vanishes 
for pairs with either (5 Q = or (3i = 0, it will in general have an important effect on the 
correlation function, especially for pairs with large \Y\. 

The effect becomes most easily apparent when we plug in some numbers and plot the 
correlation function. For simplicity, we consider a pion source with no transverse flow (v = 0) 
which freezes out instantaneously (5r = 0) with the following other source parameters: 
Rg = 3 fm, r = 4 fm/c, 8r] = 1.5, and T = 150 MeV. Restricting ourselves to pairs with 
Y = —2, Kj_ = 200 MeV, and q s = 0, we can now calculate the correlation function both 
by using the approximate analytic radii of (|1|) and by performing a numerical calculation 
directly from Eqs. (Q) and (|12"D. Comparing the results, we have found that the gaussian 
approximation of (^) with the radii of ( |T5] ) is able to describe the numerically calculated 
correlation function to within about 20% f|. Figure 1 shows a plot of the latter as a function 
of q and qi. The effect of the cross term can be seen in the form of an asymmetric ridge 
running from the peak at q Q = q\ = down to the front left where q\ > and q Q < 0. 
This kind of ridge is clearly identifiable experimentally and has in fact already been seen in 
preliminary NA35 correlation data ||17|| . 

Due to the boost invariance of the flow profile, the LCMS radii corresponding to the 
above model can be obtained simply by setting fli = and = vn t in (p!5|). Note that the 
factor Y in the second line of R 2 ol should not be set equal to zero, since it arises from the r\ 



distribution of the source ( |T2| ) which obviously breaks the boost invariance of the emission 
function. It can be verified that if a source is completely boost invariant then the cross term 
will always vanish in the LCMS frame [ l5| . However, any source with a finite size (5r] < oo) 



cannot be completely boost invariant, so it will in general feature a nonvanishing cross term 



in the LCMS, since the LCMS does not coincide with the local rest frame 18 



We would now like to suggest a better way of organizing correlation data from sources 
undergoing boost-invariant longitudinal expansion JjlJ. Returning to (|3p, let us make an 
alternative on-shell definition of the 4-vector K: 



K = (m t chY, K-t, m t shF) 



(17) 



where K 4 = ^(Pii + P2t), m t — m2 + | K t | 2 , Y = |(yi +J2), and y« is the rapidity of the 
ith particle. The resulting approximation is at least as good as the approximation we have 
been using up to now ||. This definition suggests that we express the correlation function 
in terms of q s , q D and the rapidity difference y = yi — y 2 : 



C(y, q„ q a , Y, K±) ~ 1 ± A exp[-q 2 s R 2 s - q 2 R 2 - y 2 « 2 - 2q yR 



oyj • 



(18) 



The reader should take care not to confuse the rapidity difference y with the cartesian 
coordinate y. 

The model independent expressions corresponding to (ITT]) are now given by: 



Rl = (y 2 ) 

R 2 = ([x- (K ± /m t )rch( V - Y)} 2 ) - (x - (K ± /m t )rch( V - Y)f 



2 \ r 

[m t rsh(ri — Y)\ ) — (m t Tsh(r] — Y)y 



a 



R oy = ([m t x — K±Tch(r] — F)]rsh(?7 — Y)\ — lm t x — K^rch(ri — Y))(Tsh(r) — Y)) . 



(19) 



For the model 



the radii take the much simpler form: 



R 2 S 
R 2 , 



Rl 



D 2 K l 
R* H 2 

mi 



R 



(l + (Sr])l) (5r) 2 + i(^r ( 

a 2 = m 2 t (5r ] )l[r 2 (l + u(5 V )l)+(5r) 
K ± Y 



oy 



(Sv)t (Srj)ir 2 + (5r 



(20) 



8 



The astute reader will note that the above "side" and "out" radii are identical to the LCMS 
versions of fll5f ), and that aside from a slight difference in the definition of Y, i^(LCMS) = 
a/m t and .R^(LCMS) = R oy /m t . In fact, for systems undergoing Bjorken longitudinal 
expansion, LCMS correlation functions are nothing more than approximations to fixed frame 
correlation functions in rapidity coordinates ||. 

Using the same source parameters as in Fig. 1, the effect of the cross term can be seen 
in Fig. 2 where we plot the correlator as a function of y for q D = 30 MeV. The accuracy of 
the analytic approximation (dashed) is seen by comparing it with the exact numerical result 
(solid). 

We have shown that an "out-longitudinal" (or "out-rapidity" ) cross term arises naturally 
both in a general gaussian derivation of the correlation function and in two specific gaussian 
models. Although transformed, in general the cross term persists when one switches from 
calculating momentum differences in a fixed frame to calculating them in the LCMS. Con- 
sequently, there is no reason why such a term should be excluded a priori from gaussian fits 
to experimental correlation data. Not only will the new parameter reveal more information 
about the source, its inclusion will undoubtedly increase the accuracy of the other fitted 
radii. 

ACKNOWLEDGMENTS 

We would like to thank U. Mayer, T. Csorgo, M. Gyulassy and Yu. Sinyukov for en- 
lightening discussions. Financially, this work was supported by BMBF and DFG. 



9 



REFERENCES 

[1] D. Boal, C. K. Gelbke, and B. Jennings, Rev. Mod. Phys. 62 (1990) 553. 
[2] NA35 Coll., T. Alber et al, Phys. Rev. Lett. 74, 1303 (1995). 

[3] NA35 Coll., G. Roland et al., Nucl. Phys. A566, 527c (1994); NA35 Coll., D. Ferenc 
et al., Frankfurt preprint IKF-HENPG/2-94 (1994); NA44 Coll., M. Sarabura et al., 
Nucl. Phys. A544, 125c (1992); E802 Coll., T. Abbott et al., Phys. Rev. Lett. 69, 1030 
(1992). 

[4] NA35 Coll., P. Seyboth et al., Nucl. Phys. A544, 293c (1992); NA35 Coll., D. Ferenc 
et al., Nucl. Phys. A544, 531c (1992). 

[5] NA44 Coll., H. Beker et al., CERN-PPE/94-119, submitted to Phys. Rev. Lett. 

[6] G. Bertsch, M. Gong, and M. Tohyama, Phys. Rev. C 37, 1896 (1988). 

[7] S. Pratt, T. Csorgo and J. Zimanyi, Phys. Rev. C42, 2646 (1990). 

[8] S. Chapman, P. Scotto and U. Heinz, |hep-ph/ 9409329] Heavy Ion Physics 1, 1 (1995). 



[9] S. Pratt et al., Nucl. Phys. A566 (1994) 103c; M. Lisa et al., Phys. Rev. Lett. 71, 2863 
(1993). 

[10] T. Csorgo and S. Pratt, in: Proc. of the Workshop on Relativistic Heavy Ion Physics, 
Budapest, KFKI-1991-28/A, p. 75. 

[11] G. Bertsch, P. Danielewicz, and M. Herrmann, Phys. Rev. C49, 442 (1994). 

[12] Yu. Sinyukov, in Hot Hadronic Matter: Theory and Experiment, edited by J. Letessier 
et al. (Plenum, New York, 1995), in press. 

[13] T. Csorgo, Lund U. preprint LUNFD6 (NFFL-7081) (1994), Phys. Lett. B, in press; T. 
Csorgd and B. L0rstad, Lund U. preprint LUNFD6 (NFFL-7082) (1994). 

[14] B. R. Schlei et. al, Phys. Lett. B293, 275 (1992); J. Bolz et. al., Phys. Lett. B300, 404 

10 



(1993). 

[15] A. Makhlin and Yu. Sinyukov, Z. Phys. C39, 69 (1988). 

[16] This effect arises at a higher order in (5r/)* than was studied in [F3 . 



[17] NA35 Coll., T. Alber et al, talk given at Quark Matter '95, Monterey, Jan. 6-13, 1995. 

[18] S. Chapman, J.R. Nix and U. Heinz, in preparation. 

[19] V. Averchenkov et al., Sov. J. Nucl. Phys. 46, 905 (1987). 



11 




q ou t (MeV) 

Figure 1. The numerically calculated correlation function (3) generated by the pion source (12) with 
parameters v = St = 0, Rq = 3 fm, To = 4 fm/c, Srj = 1.5, and T = 150 MeV is plotted as a 
function of qi and q for Y = —2, K±_ = 200 MeV, and q s = 0. 
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y = yi - y 2 

Figure 2. The same parameters as in Fig. 1 are used to plot the correlation as a function of 
y = yi — y2 for q = (upper curves) and q = 30 MeV (lower curves) as calculated numerically 
(solid curves) and analytically via (15) (dashed curves). 
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